Variation in displacement energies due to isospin nonconserving forces 
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For mirror nuclei with masses A = 42 — 95, the effects of isospin nonconserving nuclear forces are studied 
with nuclear shell model using the Coulomb displacement energy and triplet displacement energy as probes. It 
is shown that the characteristic behavior of the displacement energies can be well reproduced if the isovector 
and isotensor nuclear interactions with J — and T — 1 are introduced into the / 7 n shell. These forces, with 
their strengths being found consistent with the nucleon-nucleon scattering data, tend to modify nuclear binding 
energies near the N = Z line. At present, no evidence is found that these forces are needed for the upper fp- 
shell. Theoretical one- and two-proton separation energies are predicted accordingly, and locations of the proton 
drip-line are thereby suggested. 

PACS numbers: 21.10.Sf, 21.30.Fe, 21.60.Cs, 27.50.+e 



Isospin is a fundamental concept in particle and nuclear 
physics yj]. Isospin-symmetry breaking occurs in particle 
physics because of the u-d quark mass difference and the elec- 
tromagnetic effects in quarks ]2|]. In nuclear physics, nucleon- 
nucleon scattering data suggest that the neutron-neutron (nri) 
interaction is ~1% more attractive than the proton-proton 
(pp) interaction and the proton-neutron (pn) interaction is 
^2.5% stronger than the average of the nn and pp interactions 
13j-|5[]. In nuclei, the Coulomb interaction between protons 
also breaks both charge symmetry and charge independence. 
The Coulomb displacement energy (CDE), i.e. the binding- 
energy difference between mirror nuclei, is a well-known sig- 
nature of charge-symmetry breaking due to the Coulomb in- 
teraction 16J. However, it was realized |]7j] that even if the 
pairing, exchange, and electromagnetic spin-orbit terms are 
considered, the Coulomb force alone cannot account for the 
experimental CDE (known as the Nolen-Schiffer anomaly). 
There have been many attempts to resolve this discrepancy 
18[]. Shell-model calculations suggested that the isospin non- 
conserving (INC) nuclear interactions are important for under- 
standing the anomaly \%. In addition, one could also study the 
triplet displacement energy (TDE) |10y. which is regarded as 
a measure of breaking in charge independence I6J. 

The study of proton-rich nuclei is one of the frontiers in 
low-energy nuclear physics. Proton-rich nuclei with masses 
A ~ 60 — 70 are of particular interest. In this mass region, 
there are at least three so-called waiting-points along the sug- 
gested path of rapid proton capture process (the rp-process) 
1 1111 : 64 Ge, 68 Se, and 72 Kr, having equal numbers of neu- 
trons and protons (N = Z). Precise masses in the vicinity of 



the waiting-point nuclei 112141711 are required to locate the re- 
process path and to understand astronomical observations on 
the abundance of chemical elements. The concept of CDE is 
thought to be a reliable method for predictions of unknown 
masses (or nuclear bin ding energies) on the proton-rich side 
of the N = Z line fjjl, [l9D. In addition, N ~ Z nuclei with 
A ~ 80 are known 120ti to undergo dramatic changes in shape 



1 21 , |22|] with addition or removal of just one or two nucleons 
1 2311 . which would strongly influence the determination of the 
end point of the rp-process 12411 . i.e. the heaviest element that 
the rp-process nucleosynthesis may create. 
The CDE for mirror nuclei is defined as 



CDE{A, T) = BE{T, T z< ) - BE(TJ Z . 



(1) 



where T z = (N — Z)/2 is the z component of the total isospin 
T, and BE(T, T z< ) and BE(T, T z> ) are (negative) binding en- 
ergies in an isospin multiplet having the largest proton number 
(Z>) and the smallest one (Z<), respectively. For T = 1/2, the 
experimental CDE's ll 15l 12511 are shown in Fig. 02a) and com- 
pared with the Coulomb energy prediction |6J . A monotonous 
increasing trend in CDE with increasing mass number is de- 
scribed for the entire region from A = 5 to 7 1 . However, an 
overall overestimate by the calculation is seen in Fig. [Ha). 
These deviations from data can be qualitatively understood 
by the exchange effects due to the Pauli Principle, which 
keeps theprotons apart, thus weakening the Coulomb repul- 
sion la, LZD. A close examination on the curve indicates a 
zigzag behavior in these CDE's. To see the zigzag pattern 
more clearly, we introduce a quantity measuring the differ- 
ences in CDE between nuclei A and A + 2, 



ACDE(A, T) = CDE(A + 2,7) - CDE(A, T) 



(2) 



In Fig. [ljb), one clearly sees an odd-even staggering pattern. 
The Coulomb energy prediction gives only the average with 
a smooth curve. A notable exception in the pattern is seen 
for the fi / 2 -shell nuclei with masses A — 42 — 52, where the 
staggering seems to be washed out considerably. 

The TDE with T = 1 is defined with binding energies of 
triplet nuclei as 

TDE(A, T) = BE(T, T z< ) +BE(T, T z> ) - 2BE(T, T, = 0). 

(3) 
In Fig. [TJc), the known experimental TDE's 112511 are shown 
for different masses. Except for those around A = 6 and at 




FIG. 1 : Experimental Coulomb displacement energy and triplet displacement energy, (a) CDE, (b) differences in CDE between A and A + 2 
nuclei (shown as ACDE/Z), (c) TDE, and (d) differences in TDE between A and A + 4 nuclei (shown as ATDE/Z). Experimental data are taken 
from Ref. 1151 12511 . Theoretical curves from the Coulomb prediction [6|] are shown for comparison. 



A = 58, the Coulomb energy prediction disagrees strongly 
with data, particularly for those f 7 / 2 -shell nuclei where much 
enhanced TDE's are observed experimentally. We may further 
introduce a quantity measuring the differences in TDE, 



ATDE(A, T) = TDE(A +4, T) - TDE(A, T). 



(4) 



In Fig. Hid), it is seen that staggering occurs only for light nu- 
clei, but fades away for heavier ones. For nuclei starting from 
A = 30, the experimental ATDE's show a smooth behavior, 
with only A = 54 as an exception. 

Questions arise as to why in the mass region of A = 42 — 54 
the staggering magnitude in ACDE is greatly reduced, why the 
overall TDE is significantly larger than the Coulomb predic- 
tion for this mass region, and why the ATDE at A = 54 (and 
TDE at A = 58) suddenly deviates from the smooth trend. To 
find an answer, we perform state-of-the-art shell-model calcu- 
lations with inclusion of the INC interaction Hjnc in addition 
to the original isoscalar Hamiltonian Hq. For Hq, we adopt 
two modern interactions: GXPF1 A II26II with the full fp shell 
and JUN45 J27I1 with the p/5/2^9/2 model space. The total 



Hamiltonian then reads 



H = H Q +H I 



INC, 



(5) 



where Hjnc takes the form of a spherical tensor of rank two 



Hinc =H SB + Vc- 



- y v [k) 

k=\ 



(6) 



with Vc in Eq. © being the Coulomb interaction and H' 
the single-particle Hamiltonian that includes the Coulomb 
single-particle energy for protons and the single-particle en- 
ergy shifts Bis due to the electromagnetic spin-orbit interac- 
tion for both protons and neutrons with the parameters taken 
from Ref. 12811 . The Coulomb single-particle energies for pro- 
tons are taken as (all in Me V) e(0/ 7 / 2 ) =7.4, e(lp 3 / 2 ) = 7.2, 



6(0/5/2) — 7.1, and s(lpi n) = 7.3 for the fp model space, 
and e(lp3 /2 ) = 9.4, e(l/ 5/2 ) = 9.1, e(lp 1/2 ) = 10.0, and 
£(0g9/ 2 ) = 9-7 for the pfs/289/2 model space. The electro- 
magnetic spin-orbit term has been shown to play an impor- 
tant role for understanding the anomalies in the Coulomb en- 
ergy difference in 67 As/ 67 Se JH and ™Br/ 70 Se Q. The e u 
term 13 ill does not appear explicitly because this term shifts 
only the proton single-particle energies and are effectively in- 
cluded in the Coulomb single-particle energies listed above. 
V ]NC in © is the INC interaction, with k = 1 and k = 2 for 
the isovector and isotensor component, respectively. The two- 
body matrix elements with T = 1 are related to those in the 
proton-neutron formalism la LL8J through 



v*- 1 ) — V —V 

y INC ~ V PP v,h 



y (2) -v 

y INC ~ y PP 



■V m -2V, 



pn - 



(7) 



where V pp , V nn , and V p „ are, respectively, the pp, nn, and pn 
matrix elements of T = 1 . 

Calculations are performed for odd-mass nuclei with 
isospin T = 1/2, 3/2, and 5/2 and for even-mass nuclei with 
T = 1, 2, and 3, in both the fp and pf 5/289/2 model spaces. 
Because of large dimensions involved in the calculation, it is 
necessary to restrict the number of nucleons to be excited from 
the lower to the upper orbits. We have carefully checked the 
results between calculations with and without restrictions and 
found that they differ by only a few keV. 

In the GXPF1 A calculation within the full fp shell, the INC 
interaction for the fj/2 shell has terms (see Eq. d7]l) V pp = 
&»C V m = P, m V„ J =°, and Vnn = P P nVf= Q , where V J =°, 



V/ n =0 , and V/,7 are, respectively, the pp, nn, and pn pairing 
interactions for the matrix elements having a unit value. The 
parameters J3 PP = —22.5, j3,„, = 77.5, and j3 p „ = —55.0 (all 
in keV) are chosen so as to reproduce the experimental CDE 
and TDE data. Fig. [2ja) shows that the calculated CDE with 
and without the INC interaction can describe the experimen- 
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FIG. 2: (Color online) Calculations with GXPF1A and JUN45 interactions are compared with experimental data 0151 1250 for (a) CDE, (b) 
ACDE(A, T)/Z, (c) TDE, and (d) ATDE(A, T)/Z. Solid (open) symbols indicate results with (without) the INC nuclear interactions in the / 7 / 2 
shell. For comparison, results from Brown et al. 11911 and from the Coulomb prediction Jg| are also shown. 



tal data reasonably well. However, we find that only with the 
isovector and isotensor interactions can the calculation cor- 
rectly reproduce the observed reduction in staggering magni- 
tude of ACDE/Z for the mass region A = 45 — 5 1 , as shown in 
Fig. |2b), and the large experimental TDE shown in Fig. |5Jc). 
Without the INC nuclear interaction, the calculated stagger- 
ing magnitudes for ACDE/Z are clearly larger than the data 
and the TDE values are close to the Coulomb prediction, but 
smaller by about 150 keV than the experiment. In the JUN45 
calculation, the INC interaction is not included. 

The underlying physics is that inclusion of the isovector 
force in the / 7 / 2 shell modifies interactions between the nu- 
cleons. Namely, pp {nn) now becomes more attractive (less 
attractive), which results in an increase (decrease) of the pro- 
ton (neutron) pairing gap. To see its influence on ACDE di- 
rectly, we rewrite Eq. (f2} as 

ACDE(A,r)=2(-l) z >[A 7r (Z>,Z>)-A v (Z > ,Z > )], (8) 

in which A V (Z>,Z>) and A 7r (Z > ,Z > ) are the three-point odd- 
even mass differences for neutrons and protons, respectively, 
which are regarded as measures of the neutron- and proton- 
pairing gap 13211 . The occurrence of the odd-even staggering 
can then be explained by the differences between proton- and 
neutron-pairing gaps, with the factor ( — l) z> originating from 
number parity. Without the isovector force in the / 7 n shell, 
calculations give an overly strong staggering for A = 43 — 5 1 
(see Fig. |2jb)). Now with inclusion of the isovector force, an 
increasing difference between A K and A v is obtained. With 
Z> = odd (even), the factor (— l) z> in Eq. d8]l is negative 
(positive) for A = 41,45,... (A = 43,47,. . .). As compared 
to the results without the isovector force, this obviously leads 
to a decrease in ACDE for the sequence with odd Z> and an 
increase for even Z>, thus reproducing the observed reduction 
of staggering magnitudes shown in Fig. EJb). On the other 



hand, since inclusion of the isotensor force makes pn more 
attractive than the average of pp and nn, the last term in Eq. 
(0 becomes smaller, thus increasing the TDE for A = 42 — 54, 
as shown in Fig. [2jc). We note that the calculations do not 
support the apparent change in the staggering phase at A = 69 
in the experimental ACDE. This may suggest 13311 that the 
mass of 69 Br 11611 was measured for an isomer, not for the 
ground state. From the present calculations, we find that nn 
is ^0.8% more attractive than pp, and pn is ^2.5% stronger 
than the average of nn and pp. These ratios are in accord with 
those estimated from the nucleon-nucleon scattering data ]4|]. 

For the heavier mass region with A = 55 — 67, the calcu- 
lated CDE differences are also in a good agreement with the 
observed large staggering (see Fig. EJb)). The sudden drop 
in TDE at A = 58 (Fig. [2c)) an d the corresponding drop in 
ATDE at A = 54 (Fig. |2|d)) are correctly reproduced. For 
nuclei below A = 54, since nucleons occupy mainly the fan 
shell, the added INC interaction shows a significant effect, 
which correctly describes the observed large TDE, as dis- 
cussed above. For the triplet nuclei with A = 58, however, 
two nucleons occupy the p^n orbit and do not feel an INC in- 
teraction, and therefore, the TDE decreases drastically. Thus, 
in our calculation the observed sudden drop in TDE at A = 58 
may suggest that the INC nuclear interaction is less important 
for the normal-parity P3/2 and fen orbits. Differences be- 
tween the GXPF1A and JUN45 calculations are found above 
A = 69 in Fig. 2, which are attributed to the contribution from 
the g 9 n orbit. 

On the basis of the successful CDE calculation as presented 
in Fig. |2j now we try to map the proton drip-line by evaluat- 
ing one- and two-proton separation energies. According to Eq. 
|[T), we use the shell-model CDE and the observed binding en- 
ergy BE(T, T z> ) for the nucleus from the neutron-rich side to 
predict the binding energy BE(T, T z< ) for the proton-rich ana- 
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FIG. 3: (Color online) Calculated one- and two-proton separation energies for odd-mass nuclei with isospin T = 1/2,3/2,5/2 and for even- 
mass nuclei with T = 1,2,3 using the GXPF1A (left) and the JUN45 (right) interaction. In each box, the first number denotes one-proton 
separation energy and the second denotes two-proton separation energy. Thick (red) lines indicate the proton drip-line. 



logue nucleus. For binding energies above A — 82 where no 
data are available, we simply adopt the Audi-Wapstra extrap- 
olation from AME'03 112511 . The agreement between the cal- 
culated and experimental binding energies is very good within 
an rms deviation of about 100 keV. Figure [3]shows the calcu- 
lated one- and two-proton separation energies denoted in each 
box by the first and second numbers, respectively. The thick 
(and red) lines represent the proton drip-line beyond which 
the one-proton and/or two-proton separation energies become 
negative. The INC term for the fan shell is included in the 
calculation with the GXPF1A interaction. The existing data 
for 60 Ga 0, 64 > 65 As CI El, and 69 Br CI El indicate that 
these nuclei are unbound. The experimental separation ener- 
gies of 63 Ge, 67 Se, and 71 Kr IU5I1 suggest that they are bound. 
In the graph on the right, the experiments indicate that 77 Y 
and 82 Mo are bound while no evidence was found for 81 Nb 
and 85 Tc. As one can see, most of our results are consistent 
with the current experimental information. Figure |5]also sug- 
gests several candidates for proton emitters. 

In summary, we have investigated effects of the isospin 
nonconserving forces that cause characteristic shell changes 
near the N = Z line. Large-scale shell-model calculations 
were performed by employing two modern effective interac- 
tions (GXPF1A and JUN45) for the corresponding mass re- 
gions with inclusion of the Coulomb plus INC nuclear inter- 
actions. We concluded that the INC forces are important for 
the /7/2-shell nuclei, but not for the upper fp-shell. This con- 
clusion is consistent with those found in our previous papers 
II29 U3OII . No conclusion about the INC forces can currently 
be drawn for heavier nuclei with A =70 — 95. Consequently, 



we calculated one- and two-proton separation energies to map 
the proton drip-line. Our calculation provides many new pre- 
dictions for the fpg shell region up to A = 95, which may be 
relevant to the discussion of the rp-process of nucleosynthesis 
1 2411 . The results shown in the present Letter should be tested 
by future experiments on proton-rich nuclei of the heavy mass 
region. 
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